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ABSTRACT 

We  study  the  singularly  perturbed  two-point  boundary  value  problem  on  the 
interval  (-1,1) 

-eu"  +  p(x)u'  +  q (x) u  =  0,  u  (-1)  =  A,  u(l)  =  B  , 

in  which  the  coefficient  p  has  zeros  in  the  interval  (-1,1).  For  the  case 
where  p  has  precisely  one  zero  we  give  a  survey  of  our  method  of  construction 
of  an  asymptotic  approximation  of  the  solution  via  the  eigenfunction  expansion, 
and  we  show  that  "Ackerberg  -  0 'Malley  resonance"  is  identical  to  ordinary 
resonance,  namely  that  a  free  mode  in  the  solution  is  amplified  strongly  by  a 
small  divisor.  For  the  case  where  p  has  several  zeros  and  where  q  is 
identically  zero  we  construct  an  asymptotic  approximation  and  prove  its  validity 
by  a  variational  (Galerkin)  method.  The  methods  described  in  this  paper  can  be 
generalized  to  turning  point  problems  in  several  dimensions. 
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SIGNIFICANCE  AND  EXPLANATION 


\ 

V 


Consider  the  two-point  boundary  .value  problem 
r  -feu"  + 


xu 


L^-  0u  =  0,  »  u(-l)  =  A,  u(l)  =  B  , 

■  ■  rv  t  ~^C\ 

with  0  <  1.  For  certain  values  of  namely  the  eigenvalues,  the 

problem  need  not  have  a  solution;  for  ^  — >  +0  those  eigenvalues  tend  to  the 
non-negative  integers.  Although  for  values  near  such  an  eigenvalue  the 
solution  exists,  it  is  very  sensitive  to  small  changes  in  and  for  a 

long  time  it  has  been  unknown  how  to  construct  a  rigorous  asymptotic  approxima¬ 
tion  to  it.  Yet  an  asymptotic  approximation  is  of  interest  since  the  equation 
and  its  generalizations  offer  a  deterministic  model  for  the  motion  of  a  particle 
in  a  potential  field,  which  executes  a  random  walk  under  influence  of  small 
random  forces.  In  this  paper  we  construct  asymptotic  approximations  to 
solutions  of  such  problems,  in  which  the  coefficient  of  u'  has  one  or  several 

•  |X  <  r' t 

zeros  (turning  points)  and  we  prove  their  validity. 


o  1 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


THE  SINGULARLY  PERTURBED  TURNING  POINT  PROBLEM: 

A  SPECTRAL  APPROACH 

* 

Pieter  de  Groen 


1.  INTRODUCTION. 

In  a  survey  paper  [1]  on  "the  capriciousness  of  singular 
perturbations",  Wasow  gives  a  number  of  examples  of  singular 
perturbation  problems  which  display  unexpected  behaviour. 

One  of  his  examples  is  the  following  simple  looking  two-point 
boundary  value  problem  on  the  interval  (-1,1) 

eu"  -  p(x)u'  =  f  (x) ,  u(-l)  =  u(l)  =  0  (1.1) 

in  which  the  coefficient  p  changes  sign  at  x  =  0,  p(0)  =  0 
and  p(x)  0  if  x  /  0.  It  looks  natural  that  the  solution 
converges  to  a  solution  of  the  reduced  equation  pu'  =  f 
with  the  possible  exception  of  neighbourhoods  of  the  boundary 
and  of  the  point  x  =  0,  at  which  the  reduced  equation  has  a 
singularity.  Referring  to  an  unpublished  part  of  his  Ph.D.- 
thesis  (1942)  Wasow  states: 

(i)  if  p' (0)  <  0  the  solution  of  (1)  converges  on 

[-1,0)  and  on  (0,1]  to  solutions  of  the  reduced 
equation  pu'  =  f,  which  satisfy  the  boundary 
conditions  at  +1  and  -1  respectively;  the 
limit  is  discontinuous  at  x  =  0. 

(ii)  if  p' (0)  >  0  the  solution  of  (1)  diverges  in 
general. 

This  divergence  in  (ii)  is  not  "generic";  it  is  rather  excep¬ 
tional  and  is  due  to  the  lack  of  a  zeroth  order  term  in  (1.1). 
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Problems  of  type  (1.1)  have  become  famous  by  a  paper  of 
Ackerberg  and  O'Malley,  in  which  they  construct  formal 
approximations  to  the  solution  of  the  problem 

eu"  -  p(x)u'  +  q(x)u  =  0,  u(-l)  =  A,u(l)  =  B  .  (1.2) 

Their  approximations  do  not  decay  exponentially  if  p* (0)  >  0 
and  q(0)/p'(0)  is  a  non-negative  integer,  a  phenomenon  they 
have  called  "resonance".  This  paper  has  been  followed  by  a 
large  number  of  other  papers,  studying  this  subject,  mostly 
by  formal  methods. 

The  existence  of  "resonance"  for  problems  of  type  (1.2) 
can  be  proved  in  several  ways,  all  of  which  have  in  common 
the  introduction  of  an  additional  parameter,  which  performs  a 
"resolution  of  the  singularity".  Olver  [3]  constructs  an 
approximation  by  linking  together  uniform  approximations  of 
independent  solutions  of  the  equation  containing  an  addi¬ 
tional  parameter  6.  The  linking  across  the  turning  point 
yields  conditions  on  5(e)  for  "resonance"  to  occur.  Kopell 
[4]  (see  also  this  volume)  introduces  an  additional  parameter 
6  and  studies  the  continuity  with  respect  to  6  of  certain 
integral  manifolds  of  solutions,  whose  geometric  properties 
yield  criteria  for  "resonance".  Both  approaches  have  in 
common,  that  they  do  not  construct  an  approximation  of  the 
solution  of  (1.2),  but  that  they  determine  a  function  5(e) 
such  that  the  perturbed  problem 

eu"  -  p(x)u'  +  q(x)  u  +  6  (e)  u  =  0,  u(-l)  =A,u(l)  =  B  ,  (1.3) 

displays  "resonance".  In  our  approach  (5],  before  we  try  to 
solve  (1. 2),  we  first  analyse  the  associated  homogeneous 
problem 

eu"  -  p(x)u'  + q(x)u + Au  =  0,  u(-l)  =  u(l)  =  0  (1.4) 

which  can  be  recognized  as  the  eigenvalue  problem  associated 
with  (1.2).  Formally  this  looks  not  much  different  from  (1.3), 
yet  the  spectral  analysis  yields  a  quite  different  descrip¬ 
tion  of  approximations  of  the  solution  of  (1.2). 

In  this  contribution  we  shall  give  an  overview  of  the 
paper  f5)  and  we  shall  give  for  the  particular  problem  (1.1) 
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with  several  turning  points  an  analysis,  which  can  be  gener¬ 
alized  to  the  analogous  problem  in  several  dimensions.  This 
particular  problem  (1.1)  has  applications  to  exit  problems 
for  Brownian  motion  in  potential  wells,  cf.  [6] . 

2.  MOTIVATION  BY  AN  EXAMPLE. 

In  order  to  motivate  the  spectral  approach,  we  consider 
the  particular  problem 


eu"  -  xu'  +  Xu  =  0,  u(-l)  =  A,  u(l)  =  B  .  (2.1) 

Its  exact  solution  u£  can  be  expressed  in  terms  of  con¬ 
fluent  hypergeometric  functions, 


ue (x)  *  \  (A  +  B) 


F(-n;i;x2/2e) 

F(-U;i.l/2e) 


+ 


xF(j  -  jX;3/2;x2/2c) 
F (1  -  ±X;3/2;l/2e) 


(2.2) 


provided  the  denominators  are  non-zero.  Since  they  are  non¬ 
zero  for  each  fixed  X  and  e  small  enough  (non-uni formly 
depending  on  A),  we  find  the  asymptotic  formulae  for  e -*•  0, 

u  (x)  ~  A  exp{ (-1  -  x)/e>  +  B  exp{ (x  -  1) /e } ,  (2.3a) 

v» 

xf  A.  ^ 

ue(x)  ~  i  (B  +  (-l)XA)xX  (2.3b) 

+  |  (B  -  (-1) XA) exp{  -  (x  -  1) 2/2e } , 
if  X  =  0,1,2,..., 
valid  for  constant  X. 

If  one  of  the  denominators  in  (2.2)  is  zero  for  some 
X(e),  a  solution  of  (2.1)  need  not  exist,  the  homogeneous 
problem 


eu"  -  xu'  +  X(e)u  =  0,  u(-l)  =  u(l)  =  0  ,  (2.4) 


has  a  non-trivial  solution  (eigenfunction)  and  X(e)  is  an 
eigenvalue.  Actually,  it  is  well-known  that  the  denominators 
in  (2.2)  as  functions  of  X  have  denumerably  many  zeros 
XJc(e)  for  each  e  >  0;  Sturm-Liouville  theory  implies  that 
the  corresponding  set  of  eigenfunctions  e^(x,e)  is  complete. 

Let  us  now  assume  that  Z£  is  a  formal  approximation  to 
the  solution  of  (2.1),  satisfying  the  boundary  conditions 
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r 


Z  (-1)  »  A  and  Z  (1)  ■=  B  exactly,  then  we  can  insert  it  in 
equation  (2.1)  and  expand  its  residue  in  the  eigenfunctions, 

ez"  -  xZ’  «kek  . 


Sturm- Liouville  theory  then  implies 


u 


e 


a,  e, 
r  k  k 

£  *  ~  Xk(s)  * 


(2.5) 


We  see  from  this  expression  that  Z£  can  be  a  good  approxi¬ 
mation,  unless  one  of  the  denominators  in  the  right-hand  side 
is  small.  If  a  denominator  in  (2.5)  is  small  enough,  it  can 
amplify  the  corresponding  eigenfunction  in  the  expansion  of 
u£  -  Z£  so  strongly,  that  it  is  the  dominating  term  in  the 
approximation  of  u£.  Clearly,  "resonance"  in  the  sense  of 
Ackerberg  &  O'Malley  or  in  the  sense  of  Kopell  is  identical 
to  the  well-known  phenomenon  of  resonance  in  Mechanics,  namely, 
that  a  free  mode  is  strongly  amplified  by  a  small  divisor. 

Since  equation  (2.4)  is  transformed  to  Hermite’s  equa¬ 
tion  by  the  change  of  the  independent  variable  x  =  on  a 

^-interval  which  becomes  unbounded  for  e  -*•  0,  the  eigen¬ 
values  of  (2.4)  tend  to  the  non-negative  integers.  This  fact 
is  reflected  in  the  asymptotic  formula  (2.3b). 

This  example  indicates,  that,  before  trying  to  approxi¬ 
mate  the  solution  of  the  general  problem  (1.2),  we  should 
analyze  the  spectrum  of  the  associated  differential  operator 

V 

Leu  :=  -eu"  +  pu'  -  qu  (2.6) 


and  that  we  should  show  existence  of  a  solution  of  (1.2).  It 
is  this  aspect  that  distinguishes  our  approach  from  all  other 
ones. 


3.  FIRST-ORDER  APPROXIMATIONS  OF  THE  EIGENVALUES. 

For  a  study  of  eigenvalues  and  for  application  of  Sturm- 
Liouville  theory  the  operator  Le,  defined  in  (2.6),  does 
not  look  very  suitable,  since  it  is  not  selfadjoint  in  the 
usual  space  of  square  integrable  functions.  This  can  be 
amended  in  two  equivalent  ways: 


(i)  Perforin  the  Liouville  transformation, 

l  x 

v(x)  =  u(x)J£(.x),  J£(x)  :=  exp{-  !  p(t)dtl.  (3.1) 

It  transforms  the  operator  (2.6)  and  the  eigenvalue 
equation  (1.4)  to  a  self ad joint  form, 

2 

M£V  :=  -ev"  +  (Eg-  -  jP'-q)v=Xv,  v(0)=v(l)=0.  (3.2) 

(ii)  Consider  the  operator  L£  in  the  space  of  weighted 
square  integrable  functions,  in  which  the  inner 
product  is  given  by 

1  2 

(u,v)  :=  /  u(x) v(x) J  (x) dx  .  (3.3) 

-1 

With  respect  to  this  inner  product  L£  is 
selfadjoint. 

Both  methods  have  their  virtues.  We  shall  employ  the 
first  one  in  order  to  derive  a  first  approximation  to  the 
eigenvalues.  Once  the  convergence  of  the  eigenvalues  to 
well-separated  limits  is  established,  the  second  method  is 
easier  for  setting  up  an  iterative  scheme  by  which  we  can 
obtain  approximations  of  higher  order. 

The  tool  we  use  for  the  derivation  of  a  first  estimate 
is  Rayleigh's  minimax  characterization  of  the  eigenvalues  of 
a  selfadjoint  operator.  Let 

a(M£)  =  Uk(e)|k  =  0,1,...}  with  Xk  <  Xk+1  (3.4) 

be  the  set  of  eigenvalues  of  M£,  ordered  in  increasing 
sense,  then  the  k-th  eigenvalue  satisfies,  cf.  [7], 

X.(e)  =  min  max  (M  u,u)  ,  (3.5) 

K  EchJ(-1,1)  ,dimE=k+l  ueE,||u||=l  e 

where  (•,•)  is  the  usual  inner  product  in  and  ||  •  j| 

the  associated  norm. 

From  (3.5)  it  is  clear  that  each  set  of  k  +  1  indepen¬ 
dent  trial  functions  yields  an  upper  bound  for  the  k-th 
eigenvalue.  If  p(0)  is  the  only  zero  of  p  in  [-1,1], 


we  find  a  suitable  set  of  trial  functions  as  follows.  Apply 
to  equation  (3.2)  the  stretching  £  *  x  / \p' (0) |/2e  and 
expand  the  coefficients  in  powers  of  /z ,  then  we  find  the 
lowest  order  part 


d2v 

V 


+  K  V  = 


P1 (0)  + 

.2  | p *  (0)  I 


1  +  q (0) 
IP’ (0) [  j 


v  . 


(3.6) 


This  equation  has  a  decaying  solution  only  if  the  coefficient 
of  v  in  the  right-hand  side  is  a  positive  half-odd  integer, 
say  j  +  namely  (£)  exp  (-  i  £2)  .  Using  the  first  k  +  1 
of  those  as  trial  functions  for  X^  (correcting  them  such 
that  they  are  zero  at  the  boundary)  we  obtain  the  upper  bound 


-- Ve>  i 


kp'  (0)  -  q  (0)  + 

(k  +  l)p'  (0)  -  q ( 0)  +Ck£J 


if  p*  (0)  >  0  , 
if  p*  (0)  <  0  . 


(3.7) 


The  minimax  in  (3.5)  decreases  if  the  minimizing  sub¬ 
space  E  is  taken  from  the  larger  space  of  functions  which 
need  not  be  continuous  at  x  =  ±e^  and  which  have  square 
integrable  derivatives  only  in  the  subintervals  (-1,-e^), 
(-e^,+e^)  and  (e^,l).  Hence  X^  is  larger  than  the  k-th 
eigenvalue  in  the  joint  spectrum  of  the  restrictions  of  Mc 
to  those  subintervals  with  Neumann  boundary  conditions.  Since 

the  potential  term  in  M  on  the  subintervals  (-1,-e^)  and 
x  e  _  A 

(e  ,1)  xs  bounded  from  below  by  ye  for  some  y  >  0,  the 

smallest  eigenvalue  of  these  restrictions  is  larger  than 

ye  On  (-e^,e^)  the  potential  of  M  is  estimated  from 

below  by  V£, 

2 

Vx)  :=  (  Ip*  (0)  I  -  aei)2  -  i  p'  (0)  -q(0)  -  0ei  ,  (3.8) 

if  a  and  6  are  chosen  suitably.  Hence  X^  is  bounded 
from  below  by  the  k-th  eigenvalue  rrk(e)  of  the  problem 

-tu"  +  Ve(x)u  =  Xu,  u’ (±e*)  =  0  .  (3.9) 

For  this  problem  it  is  (with  the  stretching  £  =  x/e)  not 
difficult  to  show,  that  its  eigenvalues  satisfy 
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f 


(3.10) 


kp'(0)  -  q  (0)  +  Ofe*)  if  p' (0)  >0 

Ve)  =  h 

K  (k  +  l)p'(0)  —  q  ( 0 )  +0(ec)  if  p' (0)  <  0  . 

So  we  have  given  a  sketch  of  the  proof  of  the  following 
theorem: 

Theorem  1:  The  eigenvalues  of  L£  satisfy  for  e  ■+  0  the 
estimates: 

kp'(0)  -  q (0)  +  0(ei)  if  p'(0)  >0, 

\(t)  =  ,  (3.11) 

K  (k  +  l)p'(0)-q(0)+0(e£)  if  p'(0)<0. 

For  details  of  the  proof  we  refer  to  [5]. 

If  p  has  several  simple  zeros  in  the  interior  of  the 
interval,  every  zero  produces  a  set  of  eigenvalues,  which 
satisfies  the  analogue  of  (3.11).  If  p  has  a  zero  at  the 
boundary,  this  zero  produces  "half"  of  the  above  set  of  eigen¬ 
values,  namely  only  the  ones  with  odd  index,  the  even  ones 
being  ruled  out  by  the  boundary  condition.  The  proof  is 
analogous  to  the  one  above.  For  every  zero  of  p  we  now 
construct  a  set  of  He rmite- functions  as  trial  functions;  this 
yields  an  upper  bound.  For  a  lower  bound  we  make  the  same 
type  of  restrictions  to  0 (e^) -neighbourhoods  of  the  zeros  of 
p  as  above. 

4.  APPROXIMATIONS  OF  HIGHER  ORDER  TO  THE  EIGENVALUES  AND 
EIGENFUNCTIONS. 

In  order  to  obtain  approximations  of  higher  order,  we 

return  to  the  original  operator  L£,  which  is  selfad joint 

with  respect  to  the  weighted  inner  product  (3.3).  We  assume 

now  that  p' (0)  is  positive,  such  that  the  weight  function 
2 

J£  takes  its  maximum  at  x  =  0  and  is  exponentially  small 
(relative  to  the  maximum)  outside  a  small  neighbourhood  of 
this  point.  This  implies  that  our  norm  hardly  notices  errors 
of  an  approximation  if  they  are  only  of  polynomial  order  in 
e.  Therefore  the  algorithm  we  devise  here  for  the  approxima¬ 
tion  of  the  eigenfunction  does  not  care  for  errors  outside  a 
neighbourhood  of  the  point  x  =  0. 

For  the  construction  of  a  formal  approximation  of  the 
k-th  eigenvalue  A^(e)  and  associated  eigenfunction 
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ek(x,e)  we  perform  the  stretching  £  :=  x/p'  (0)/2e  in  th 
eigenvalue  equation  (1.4),  we  insert  the  formal  expansions 


2p(C/2e/p'  (0) )//2ep'  (0)  =  2C  +  £  PHe?  C 

i=l  1 


Ji,.i+1 


(4.1) 


2q(C/2e/p'  (0)  )/p*  (0)  =  l  q.e*V 

i=0  1 


(4.2) 


Xk(e)  =  kp*  (0)  -  q ( 0)  +  |  p*  (0)  Xkje; 


(4.3) 


.(C/2e/p*  (Q).e)  =  £  ekj  (?) 


.43 


3-o 


(4.4) 


and  we  collect  equal  powers  of  e.  This  results  in  the 
recursive  system  of  equations  (e'  =  de/d£) 

ekj  "  2?ekj  +  2kekj 

i=l  (Pi5  ^  +  q +  Xki)ek,j-i-l  ' 


(4.5) 


from  which  the  unknown  functions  e 
are  determined.  Clearly,  e. 


kj 


and  the  coefficients 
satisfies  Hermite's 


ki  — — —  — -*jr»  -k0 

equation;  all  its  solutions  are  exponentially  growing  except 


the  Hermite  polynomial  Hk-  Hence  we  may  choose  evn  =  H 


kO 


V 


The  next  term  ekl  satisfies  Hermite's  equation  with  a 
polynomial  right-hand  side. 


e 


kl 


-  2?e^  +  2ke 


kl 


*  Pi«2«i  + 


(qLS  + 


xki>Hk 


(4.6) 


Again  its  solution  is  exponentially  increasing,  unless  it  is 
a  polynomial.  The  particular  equation 

y"  -  2£y'  +  2ky  =  H ^  (4.7) 

has  the  polynomial  solution  y  =  H^/(2k  -  2 j ) ,  provided 
k  ^  j.  If  k  =  j  every  solution  of  it  is  exponentially 
increasing.  Writing  the  right-hand  side  as  a  sum  of  Hermite 
polynomials,  we  see  that  a  unique  coefficient  Xk^  exists, 
such  that  this  right-hand  side  does  not  contain  a  multiple  of 
Hk*  More  generally,  in  every  step  of  the  recursion  (4.5)  we 


can  determine  a  unique  which  suppresses  the  presence  of 

exponentially  growing  terms  in  the  solution.  Finally  we  show 

that  the  odd  coefficients  A^  2i-l  ^  =  1»2...)  must  be  zero. 

If  k  is  even  (odd) ,  then  is  an  even  (odd)  function  and 

the  right-hand  side  of  (4.6)  is  odd  (even)  provided  A^  =  0 

and  cannot  contain  a  non-zero  multiple  of  H^.  Analogously, 

if  j  is  odd  and  A^  =  0  (i  =  l,3,...,j),  then  the  right- 

hand  side  of  (4.5)  is  odd  (even)  is  k  is  even  (odd). 

In  order  to  prove  validity  of  those  expansions,  we 

define  the  partial  sums  A.  .  and  E.  .  by 

K  J  K3 


j-1 


Akj(e)  :=  kp*  (0)  -  q(0)  +  |  p»  (0)  l  A 


2  j-1 


L  k  2i 
i=l  K,<ix 


Ek;.(x,e)  :=  p(x)  .  I  ev.  (x/2e/p '  (0) )  e*1  , 


(4.8) 


i=0 


ki 


where  p  is  a  C00  cut-off  function, 


p(x)  =1  if  |x|  <  j,  p(x)  -  0  if  |x|  >  |  .  (4.9) 


The  above  construction  implies 

II'1  -  Vex3II„=  0(eillEkjllw>  •  «-10> 

Expanding  E,  .  in  the  true  eigenfunctions  and  using  the 

KJ 

initial  estimate  (3.11)  we  easily  find  (cf.  [5]): 

Theorem  2:  The  eigenvalues  of  L£  admit  the  asymptotic 
expansion 


1  j;1  i  -j 

Ak(e)  =  kp’  (0)  -  q  ( 0)  +|P’  (0)  J  \f2ie  +  0(e3)  ,  (4.11) 


(e  ■*»  0) ,  Vj  , 


and  the  associated  eigenfunctions  satisfy: 

HEkj  -  "kAllw-  °<£3llEkjll„>'  «-12> 

where  “k3  !=  Hek  ^ HEkJ  »w  • 

Obviously,  the  estimate  (4.12)  in  the  weighted  norm 
can  give  good  pointwise  estimates  of  the  error  only  in  a 
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neighbourhood  of  diameter  0{/c)  of  the  maximum  of  the 
weight  function.  Such  a  pointwise  estimate  can  be  obtained 
from  the  following  variant  of  Sobolev's  inequality, 

u2(x)J2(x)  =  J  ^  (u2(x)  J2(x))dx  =  (4.13) 

a 

*  2<u,u')w  -  (Pu,u)w/e  <  C||u||2/e  +  e||u'  ||2  < 

i  c1llull^/£  +  I! u II w !! Leu !l Vu  e  hJ(-i,d 

where  we  used  the  identity 

||u'||2=  (-u"  +  pu'/e,u)w  =  (L£u  +  qu^^e  .  (4.14) 


These  estimates  (4.12-13)  imply  the  pointwise  estimate 

jEkj  (x,e)  -  ukjek(x,e)  |  Vx  e  (-y/e,y/c)  .  (4.15) 


The  approximations  Ek^  to  the  eigenfunction  ek  we 
have  constructed  up  to  this  point,  are  in  fact  the  internal 
boundary  layer  terms,  valid  only  in  a  neighbourhood  of  the 
turning  point  x  =  0.  In  order  to  construct  a  uniformly 
valid  approximation  we  have  to  match  the  internal  layer  terms 
to  regular  expansions  valid  in  (-1,0)  and  (0,1).  The 
integration  constants  are  uniquely  determined  by  the  matching 
(cf.  (51).  These  regular  expansions  are  matched  to  the 


boundary  conditions  in  ordinary  boundary  layers.  If  we  start 

with  the  normalization  e^  =  Hk,  the  approximation  is  of 

order  unity  in  an  0 (/e) -neighbourhood  of  x  =  0  and  hence 

-k/2 

of  the  order  0(e  '  )  globally  on  the  interval. 

The  validity  of  this  formal  approximation  can  be  proved 
by  common  barrier  function  techniques.  Since  we  already  have 
a  good  approximation  in  an  0 (/?) -neighbourhood  of  x  =  0, 
we  can  restrict  the  problem  to  subintervals  (-l,-Y/6f)  and 
(Y/e,l)  for  a  suitably  chosen  y  >  0,  where  the  lowest 
order  term  of  the  regular  expansion  multiplied  by  log(x)  is 
a  good  barrier  function,  cf.  [5  &  7]. 
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L 


Remarks : 

(i)  If  p’  (0)  <  0  we  can  prove  the  analogue  of 

2  * 

Theorem  2  by  considering  the  L  -adjoint  L£  of  L£, 
L£u  i-  -eu"  -  (pu)  '  -  qu,  u  c  hJ  .  (4.16) 

It  has  the  same  eigenvalues  as  Le  has,  and  the 

sign  of  p  is  reversed,  such  that  we  can  apply 

* 

Theorem  2  directly.  The  eigenfunctions  e,  of 
*  K 
L£  and  efc  of  L£  satisfy  the  relation 

ekJe2  =  ek  (4*17) 

(ii)  In  general  there  is  no  reason  why  the  asymptotic 
series  for  X^fe)  should  converge,  hence  it  is 
impossible  to  determine  approximations  of  it  with 
exponentially  small  errors,  unless  p  and  q  are 
analytic. 

(iii)  If  p  has  several  zeros,  the  above  method  for 

determining  asymptotic  power  series  for  the  eigen¬ 
values  works  equally  well.  It  is  not  yet  clear 
how  we  can  construct  uniformly  valid  approximations 
to  the  eigenfunctions  on  the  whole  interval. 

5.  EXPONENTIAL  DECAY  AND  RESONANCE  (THE  CASE  p' (0)  >  0) . 

If  no  eigenvalue  of  L£  is  zero,  a  solution  for  problem 
(1.2)  exist.  For  the  case  p' (0)  >  0  we  shall  construct  an 
approximation  following  the  suggestion  of  the  example  in 
Section  2. 

For  problems  of  this  type  boundary  layers  are  generally 
expected  at  both  end  points  of  the  interval,  since  the  singu¬ 
lar  solution  of  L£u  =0  is  increasing  near  x  =  +1  and 
decreasing  at  x  *  -1.  Hence  we  can  by  usual  methods  con¬ 
struct  a  foimal  approximation  Z£,  which  consists  of  boundary 
layer  terms  at  x  =  ±1  only,  which  is  exponentially  small  in 
the  interior  of  the  interval  and  satisfies  for  some  k  the 


f 


and  the  boundary  conditions 


2  (-1 )  «  A,  2  (1)  =  B  . 

v  t* 


(5.2) 


We  easily  find  the  lowest  order  term: 

2£(x)  =  A  exp (p (1)  (x  -  l)/£)  (5.3) 

+  B  exp(p(-l)  (-1  -  x)/e)  +  0(e)  . 

Following  the  suggestion  of  formula  (2.5)  we  expand  the 
residue  in  the  true  eigenfunctions, 


LeZe  =  E  6iei' 
e  e  j=0  3  3 


(LeZe'ellw 


Bj(e)  :  (e.,e^) 


3  w 


Hence,  if  A^fe)  -*•  0  (e  ->  0)  , 


IK  -  \  -  skvxkllw=  “(IKMIJ  • 


e  e  "  w 


(5.5)* 


Since  the  weight  function  in  the  weighted  inner  product  is 
exponentially  snail  outside  a  neighbourhood  of  x  =  0,  the 
coefficient  3^  is  exponentially  small  (in  [5]  we  compute  it 
with  a  relative  error  of  order  0(/e)).  If  the  asymptotic 
series  (4.3  &  8)  contains  non-zero  terms,  then  obviously  the 
eigenfunction  component  in  (5.5)  is  exponentially  small 
and  its  presence  in  the  approximation  is  not  noticed.  On  the 
other  hand,  if  all  terms  of  the  asymptotic  series  for  Ak(e) 
are  zero,  and  if  the  coefficients  p  and  q  are  of  class 
C°°  only,  we  cannot  determine  by  asymptotic  expansions  in 
powers  of  e  whether  Ak(e)  7*  0  for  e  ^  0,  nor  can  we 
decide  whether  8^ (e) /A^ (e)  converges  to  a  definite  limit 
(if  the  denominator  is  non- zero) .  Hence,  for  general  C™ 
coefficients  we  are  unable  to  determine,  whether  a  solution 
exist  for  each  small  enough  z  and  whether  it  converges  for 
e  -  0,  if  | Ak(e) |  <  eN,  VN. 

If  we  have  analytic  coefficients  p  and  q  and  if  the 
asymptotic  series  for  ^(e)  vanishes,  then  we  can  apply 
uniform  reduction  theory  and  connection  formulae.  If  the 
real  interval  (-1,1)  is  contained  in  the  smallest  of  the 
disks  in  complex  plane,  in  which  the  power  series  expansions 
of  p,  x/p  and  q  are  convergent,  then  it  can  be  shown. 


cf.  Sibuya  [9],  that  the  equation  L£u  *  0  has  a  solution 
which  converges  uniformly  on  [-1,1]  to  a  non-trivial  solu¬ 
tion  of  the  reduced  equation  pu*  -  qu  =  0;  i.e.  it  can  be 
shown  that  the  equation  L£u  =  0  shows  "resonance"  in  the 
sense  of  Kopell  [4] ,  in  that  case.  From  this  particular 
solution  plus  a  boundary  layer  approximation  we  can  construct 
easily  an  approximation  to  the  solution  of  (1.2).  By  adding 
boundary  layers  at  both  sides,  we  can  also  find  a  good  esti¬ 
mate  of  the  k-th  eigenvalue,  which  is  of  the  same  order  as 

is.  We  remark  that  it  is  not  known,  how  to  find  a  satis¬ 
factory  estimate  of  A^ts)  if  the  condition  about  the  radii 
of  convergence  is  not  satisfied. 

6.  A  MULTIPLE  TURNING  POINT  PROBLEM. 

6. a)  Position  of  the  problem.  We  now  return  to  the  example 

(1.1)  of  Wasow;  however,  instead  of  an  inhomogeneous  equation 
we  study  the  case  with  inhomogeneous  boundary  conditions, 

L£u  :=  -eu"  +  p(x)u'  =  0,  u(-l)  =  A,  u(l)  *  B  ,  (6.1) 

in  which  the  coefficient  p  has  several  zeros  in  the  interior 
of  the  interval  (-1,1).  This  problem  can  easily  be  solved 
exactly: 

u£(x)  -  A  +  (B  -  A)7(x)A  (1) ,  (6.2) 

x  t 

y(x)  ?=  /  exp {/  p(e)ds/e}dt 

-1  -1 

From  this  formula  it  is  easily  seen  that  the  solution  has 
transition  layers  at  the  absolute  maxima  of  a  primitive  of  p 
and  is  almost  constant  elsewhere. 

Nevertheless,  the  study  of  this  problem  is  interesting 
from  the  point  of  view  of  its  analogue  in  several  dimensions, 
where  an  exact  solution  does  not  exist.  In  [6]  and  [10]  a 
variational  method  is  described  which  yields  a  formal  approx¬ 
imation,  both  in  one  and  in  several  dimensions.  In  this  sec¬ 
tion  we  shall  sketch  a  proof  of  the  validity  of  this  varia¬ 
tional  method  in  one  dimension,  without  using  the  information 

(6.2)  we  have  about  the  exact  solution.  The  proof  carries 
over  to  several  dimensions,  but  is  more  complicated  there. 


For  a  simplification  of  the  presentation  we  shall  make 
the  following  assumptions 

(i)  all  zeros  of  p  are  simple,  p(-l)  >  0  and  p ( 1 )  <  0, 

(ii)  p  has  2n  +  1  zeros;  the  zeros  at  which  p'  is 

negative  are  denoted  by  i  =  0,...,n,  and  the 

zeros  at  which  p*  is  positive  are  denoted  by 
a^,  i  =  l,...,n;  moreover  we  define  :=  -1 
and  <*n+i  *=  1/  such  that  we  have 

°0  <  30  <  °1  <  81  •**  <  0n  <  an+l  ‘  (6,3) 

(iii)  A  »  0  and  B  =  1. 

There  is  nothing  deep  in  those  assumptions,  they  relieve 
vis  of  having  to  write  down  several  alternatives  in  most  of 
the  formulae.  E.g.  a  boundary  layer  at  the  boundary  points 
±1  in  general  has  a  different  structure  than  an  internal 
layer  has;  the  assumption  ±p(±l)  <  0  precludes  layers 
at  *1. 

6.b)  The  variational  formulation.  The  variational  approach 
of  [6]  and  [10]  consists  essentially  of  the  following.  The 
operator  L£  is  sel fad joint  with  respect  to  the  weighted 
inner  product  ( • ,  • )  , 

1 

(u,v)  :«  /  u(x)  v(x)exp{-P(x)/e}dx  ,  (6.4) 

w  -1 

x 

p(x)  :=  /  p  (t)  dt 
*1 

where  the  index  I  is  chosen  such  that  P(x)  <_  0  for  all 
x  e  [-1,1],  Equivalent  to  problem  (6.1)  is  the  variational 
problem  to  find  u  e  Hg(-l,l)  which  satisfies  the  boundary 
conditions  u(-l)  =  A  and  u(l)  =  B  and  the  variational 
form 

Bg(u,v)  :«  e(u',v')w  “0,  V  v  e  Hg(-l,l)  .  (6.5) 

When  we  try  to  find  approximations  to  the  solution  of  (6.5), 
we  see  at  once  that  the  best  candidates  for  approximations 
are  among  those  ones  which  are  almost  constant  everywhere. 


except  possibly  near  the  minima  0^(i  *  0,...,n)  of  the 
weight  function,  where  internal  transition  layers  may  arise. 
In  [6]  and  [10]  a  formal  approximation  is  constructed,  which 
is  the  solution  of  the  restriction  of  B_  to  a  finite  dimen- 
sional  subspace  consisting  of  such  functions.  We  shall  prove 
here  the  validity  of  such  a  method  by  an  argument  coitmon  for 
Galerkin  methods.  For  it  we  have  to  find  a  subspace  which 
contains  good  approximations  of  the  solution  and  of  the  test 
functions;  moreover,  we  have  to  find  a  suitable  positive 
lower  bound  for  the  bilinear  form  Be. 

6.c)  Construction  of  a  formal  approximation.  In  order  to 

construct  the  subspace  in  which  we  look  for  an  approximation 

of  the  solution  u  of  (6.1)  we  first  construct  error 

e 

function-like  formal  approximations  to  the  expected  internal 
transition  layer  terms  at  0 .  (i  =*  0,  .  . . ,  n)  and,  thereafter,  we 
show  that  a  linear  combination  thereof  may  yield  a  good 
approximation  of  u£.  To  this  aim  we  define  the  functions 
pi»*i>  <i»i  and  xLi 

is  a  C°°  cut-off  function  which  satisfies 


P^x) 


if  lx-0.1<|h, 
if  |x  -  ei|  >|h. 


h  :  =  min  0.  -  a .  ; 
i  1  1 


(6.6) 


<fi .  is  an  approximation  to  the  non-constant  part  of  P  near 
Bi' 

*A(x)  :=  l  Pm  (f^)  (x  -  e.)Vj!  -  «i(x  -  0-T*  (6.7) 

j=2 

2k 

in  which  is  chosen  such  that  v ^  (x)  £  P(x)  -  P  (6^)  -  (x  -  0^)  , 

Vxt  is  an  approximation  of  the  error  function- like 
transition  layer  at  0^, 

x 

^(x)  :=  yi(e)  /  (t) exp^  (t)/e }dt  ,  (6.8) 

•00 

in  which  y^e)  «  (-P*  (8i)/2ire)  ^  (1  +  0(e))  is  such  that 
♦A(x)  »  0  if  x  £  and  ^(x)  -  1  if  x  ai+i*  is 

the  function 


XA(x)  :»  “  ^(x),  (i  -  l,...,n)  (6.9) 
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which  is  nearly  the  characteristic  function  of 
employed  in  [6]  as  test  function  in  the  variational  form. 
The  construction  implies  that  L£i^  is  zero  outside 
(cu_^#ou)  and  that  it  satisfies  the  estimates 


I Le (x )  I  -  wiCix  ’  Bi|2k“1exp{*i(x)/e} 
<  ei“1/2k  c  exp{  (P(x)  -  P^M/e}  , 


(6.10a) 


(6.10b) 


inside  the  interval. 

A  linear  combination  of  these  functions  is  to  be  an 
approximation  of  the  solution  of  (6.1).  This  linear  combina¬ 
tion  is  determined  by  the  discretization  of  the  variational 
form  (6.5).  As  test  functions  we  take  the  subspace  E, 

E  spanfx^  .  •  •  ,Xn> 

and  we  seek  a  solution  v„  in  the  linear  manifold  ip  +  E, 

e  rn 

which  consists  of  all  linear  combinations  of  \p^  (i  =  0, . . .  ,n) 
that  can  satisfy  the  boundary  conditions  u(-l)  =  0  and 
u(l)  *  1, 


MM  AM 

v_  :*  l  with  J  r\.  ■  1  • 

e  i-0  1  1  i«0  1 


(6.11) 


This  solution  of  the  discretized  variational  form  has  to 
satisfy 


WX>  =  °'  v  X  «  E  . 

Taking  a  basis  in  E  we  find  the  set  of  n  equations 


°  *  w*!1  *  * 


(6.12) 


-  "i-iiin-iii;  -  "iiiniiw 

From  these  equations  and  the  side  condition  £  =  1*  the 

coefficients  can  be  solved;  it  is  easily  seen  that  they  are 
all  positive.  Approximate  evaluation  of  the  integrals  yields 
the  equations 

n^iU^expt-Pta,  ,)/e)  «  n4y,exp(-P(B,  )/e)  (1  +  0(e) )  . 


(6.14) 


Since  P($i)  <0  and  P(8j)  =  0,  this  Implies 

n ±  »  OtexpfPlB^/e) ) ,  i  =  0,...,n  . 

6.d)  A  lower  bound  for  the  bilinear  form.  A  lower  bound  for 
B£  is  derived  with  the  aid  of  the  maximum  principle  and  a 
suitable  barrier  function.  We  shall  employ  the  following 
generalization  of  [8,  Theorems  14  and  17] : 

Lemma  1:  If  a  positive  continuous  function  W  exists,  which 
2 

is  C  except  in  the  points  {x^,...,xn)c  (-1,1)  and  which 
satisfies 


L£W(x)  >  AW (x) 

V  x  *  xi 

(i  ■  1, . . . ,n)  , 

(6.15a) 

W'  (xi  -  0)  >  W* 

(xt  +  0) 

(i  *  l,...,n)  , 

(6.15b) 

then  the  smallest  eigenvalue  of  L  is  not  smaller  than  A. 

2  e 
If  moreover,  u  e  C  satisfies 

u(±l)  =  0  &  |L  u(x)  |  <  L  W(x)  ,  V  x  ?  x.  ,  (6.16) 

then  u  is  bounded  by  W, 

|u(x)  |  <  W(x),  V  x  e  (-1,1]  . 

Proof:  Assume  that  u  -  W  has  a  positive  maximun  at  an 
interior  point  x  =  a  /  x . (i  =  l,...,n),  then  L~(u  -  W)  is 

X 

positive  at  a,  which  contradicts  (6.16).  If  u(x^)  -W(x^) 
is  a  local  maximum,  then 

u'  (x^)  -  W’  (x^  -  0)  >  0  and  u'  (x^)  -  W  (x^  -  0)  <  0  ; 

this  contradicts  (6.15b).  Hence  u  -  W  has  no  interior 
maximun  and  is  negative  at  the  boundary,  and  thus  negative 
everywhere.  Likewise  -u  -  W  is  negative  everywhere.  For 
the  assertion  on  the  eigenvalue,  we  assume  L£u  =  Xu  with 
X  <  A;  from  L£ (u/w)  <0  it  then  follows  analogously  that 
u/w  has  no  interior  maxima  nor  minima  and  hence  that  it  is 
zero  everywhere.  O 

A  suitable  barrier  function  is  the  function  W, 


I 


(6.17) 


W(x)  »  minj/  (1  +  t)  exp{P(t)/e)dt,  (6.) 

x  ^ 

/  (1  -  t) exp{P (t) /e }dtj  ; 

it  satisfies  the  correct  jump  condition  across  the  discon¬ 
tinuity  of  W'  and  outside  this  point  it  satisfies 

L£W(x)  =  e  exp{P (x)/e)  . 

A£  is  the  minimum  of  this  residue  divided  by  W.  Let  l 
the  index  for  which  the  expression 


P(a.)  -  min{max  P(8.)»  max  P($.)> 
1  j<i  1  j>i  1 


(6.18) 


takes  its  minimal  value,  and  let  m  and  r  be  the  indices 
at  which  {P(0^)|i  <  4}  and  {P(B^)|i  >  Jt)  respectively 
take  their  maximal  value.  Obviously,  r  =  I,  since  P  has 
an  absolute  maximvm  at  ST.  Standard  computations  show: 


A  *  min  e  exp(P  (x)/e)/W(x)  = 
x 


(6.19) 


=  (-ep*  (Bm)/2Tr)Jexp{P(ctjl)/e-P(Bm)/eHl  + 0(e))  . 


It  is  easily  seen  from  the  minimax  criterion  (3.5)  and  the 

trial  function  il>  -  ij)  that  Ce  ^A_  is  an  upper  estimate 

m  r  £ 


for  the  smallest  eigenvalue.  This  shows  that  A£  is  a 
nearly  optimal  lower  bound  for  B£. 

6.e)  Approximation  properties  of  the  subspace.  The  second 


point  on  which  the  proof  of  the  validity  of  the  variational 
approximation  v£  is  based,  is  -  as  usual  in  Galerkin 
methods  -  that  the  approximate  solution  space  ip  +  E  con¬ 
tains  a  satisfactory  approximation  of  the  true  solution  u£ . 
We  shall  show  that  such  an  approximation  is  given  by  Ue, 


U£  ^n  +  E  ue(cti)xi  =  E  ^i 
e  n  ±ml  £  i  i  i  1 


(6.20) 


where 


ue^ai+l^  “  uc^ai^*  Since  it  follows  from  the 


maximun  principle  that  u£  is  monotonely  increasing,  all 
are  positive  and  smaller  than  or  equal  to  one.  Their 
size  can  be  estimated  much  better. 


I 


a 


Let  us  during  this  paragraph  assume  that  P  has  only 
one  absolute  maximum  in  the  interval  (-1,1),  viz.  Bj. 
Define  the  barrier  function  WQ, 


WQ(x) 


x 

/  (8t  -  £)exp(P(5)/e)d£/e,  if  x  <  BT  , 

-1  1  “  1 

(6.21) 

1 

Y  /  (S  -  BjJexptPtU/eJdC/e,  if  x  >  8  , 


where  y  =  (1  +  0(e))  is  chosen  such  that  Wq  is  continuous 
at  Bj.  As  in  (6.17)  this  function  satisfies 

Lew0  =  expfPOO/e}  ;  (6.22) 


moreover,  its  derivative  is  continuous  at  Bj  and  it  has  a 
maximum  there  of  order  unity.  More  generally  we  have 

W_(a.)  =  0(max  e“^exp(P(Bj)/e) ) ,  if  i  <  I  ,  (6.23) 

P  1  j<i  3 


and  an  analogous  formula,  where  maximum  is  taken  over  j  i, 
if  i  >  I.  Let  now  Y  be  a  linear  combination  of  \p^. 


n 


*  -  I  Vi' 

i=0  1 


I.  "  1 


i=0 


(6.24a) 


in  which  the  orders  of  the  coefficients  are  prescribed  by  the 
conditions 


*  0(exp(P(Bi)/e) ) ,  i  =  0,...,n  .  (6.24b) 

According  to  formula  (6.10b),  this  linear  combination  satis¬ 
fies  the  estimate 


L£Y  <  Ce*"1/2kexp{P(x)/e}  . 

Since  u  -  Y  has  zero  boundary  values,  we  can  apply  the 

e  4-l/2k 

second  half  of  Lenina  1  with  barrier  function  e  WQ. 

This  implies  that  u£  -  Y  is  bounded  by  e”^2kWg,  and  in 

particular  this  shows 

h  ’  VW  -  W  -  o«r1/2kw0<ai)>  . 
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At  at  least  one  point  at  the  left-hand  side  of  6j  we 

thus  obtain  an  estimate  for  of  the  form 

q  -  0(e"1/2kexp(P(8i))  .  (6.26 

and  analogously  at  the  right-hand  side  of  8j. 

This  process  sketched  above  we  can  apply  to  each  sub¬ 
interval  (0^,0^)  of  (-1,1)/  which  contains  only  one 
absolute  maximum  of  P.  Hence,  starting  with  intervals  con¬ 
taining  the  largest  maxima,  we  can  go  down  step  by  step,  in 
each  step  establishing  an  estimate  for  at  least  two  coeffi¬ 
cients  and  loosing  in  each  step  a  factor  of  order 

— l/2k  1 

e  '  .  Starting  with  a  sufficiently  large  k,  we  obtain 

in  this  way  the  estimate 

■  0(e“^exp(P(8i) ) ,  i  *  0,...,n  .  (6.27) 

For  an  estimate  of  the  global  error  of  the  approximation 

U„ ,  which  we  have  constructed  in  4>  +  E,  we  consider  a 

e  n 

subinterval  (cu,ai+1)  .  At  both  endpoints  UG  -  u£  is  zero 
by  definition;  moreover,  the  estimates  (6.10a)  and  (6.27) 
imply 

||Le(U£  -  u£)||Wfi<  Cek_2  ,  (6.28a) 

where  II  •  II  .  .  is  the  restriction  of  the  norm  to  the  sub- 
interval  Since  this  subinterval  does  not  contain 

in  its  interior  any  zero  of  p  at  which  p'  is  positive, 
the  smallest  eigenvalue  of  the  restriction  of  L£  to  this 
subinterval  (with  Dirichlet  boundary  conditions)  is  bounded 
away  from  zero.  Hence  (6.28a)  implies  that  (U£  -  u£)  is 
of  the  same  order.  Adding  up  over  all  subintervals,  we  thus 
find: 

llUe  -  ucHw<  Cek~2  .  (6.28b) 

6.f)  Error  bounds  for  the  variational  approximation.  Having 
now  a  lower  bound  for  B£  and  a  good  approximation  in  <J^  +  E, 
we  can  apply  the  usual  Galerkin-arg  ument.  From  (6.5)  and 
(6.12)  we  find 

B£(ue  -  ve,x)  *  0,  V  x  «  E  . 
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Adding  and  subtracting  U£  we  find 

Be(Ue  -  VE,X>  -  B£(UE  -  uE,X)  . 

Since  U£  -  v£  is  an  element  of  E,  we  may  choose  x  equal 
to  this.  Hence  we  find 

Vue  *  V°e  -  V  *  VU£  '  Vuc  '  V  * 

"  2Be<U£  -  Vut  '  V  '  V°e  -  v£'ue  *  V  i 

<  Be(De  -  Vue  -  u£)  , 


where  we  used  the  estimate 

B£  (u,v)  <  (B£(u,u)B£  (v,  v) )  *  <  |-Be  (u,u)  +  ^-Be  (v,  v)  . 

Together  with  the  lower  bound  (6.19)  this  implies 

UUE  -  vj|2<  Ce2k"9/2exp{|  P(Bm)  -  |  Pte*)}  . 

This  does  not  at  all  look  like  a  satisfactory  estimate.  How¬ 
ever,  if  we  evaluate  |]u  -  v  ||  approximately,  we  find 

C,  fc»  W 

Hu  -  v  W2  = 
m  e  enw 

n  i-1  -  i 

•  l  l  (C.,  -  tk)  ^  (2ire/p’  (a. ) ) *exp{-P (a. )/e}  (1+  0(e)), 
1*1  j*0  3  3  1  1 

hence  at  least  at  we  find  a  good  estimate,  namely 

IW  '  ve  (°t>  1 2  ’  £  -  lj>2 

<  Ce2k~9/2exp{P(Bm)/c}  . 

This  estimate  we  get  precisely  at  the  point  in  the  interval, 
that  "generates"  the  smallest  eigenvalue  of  L£.  If  there 
are  several  of  such  points,  i.e.  if  (6.18)  does  not  determine 
a  unique  i,  then  at  all  those  points  such  an  estimate  holds. 

Now  we  can  split  the  interval  in  (at  least)  two  sub¬ 
intervals,  (-l,o4)  and  (aA,l)  and  restrict  the  problem 
(6.1)  to  both  subintervals.  Since  is  now  a  boundary 

point  for  both  problems,  it  does  not  generate  an  eigenvalue 
that  tends  to  zero,  see  Theorem  1  and  the  comments  that 


#  - 


1 


follow  it.  Hence,  applying  the  same  proof  as  before,  we  get 
error  estimates  at  two  other  zeros  of  p,  which  in  the 
restricted  problems  generate  the  smallest  eigenvalues.  So 
we  can  go  on  until  we  have  obtained  error  estimates  of  type 
(6.29)  for  all  unknowns  in  the  variational  approximation. 
Thereafter  we  can  easily  derive  error  estimates  for  all  other 
points  of  the  interval  by  the  maximum  principle.  So  we  have 
finally  proved: 

Theorem  3:  The  approximation  v£,  generated  by  the  dis¬ 
cretized  variational  form  (6.12)  satisfies  the  error  estimate 

|u£(x)  -  v£(x) |  <  Cek‘9/4  .  (6.30) 

Remarks: 

ij  If  p(-l)  <0,  we  have  to  add  to  the  set  of  trial 
functions  {iJk)  a  trial  function  whose  exponential  part  is 
of  the  form  exp(p(-l)  (x  +  l)/e)  and  which  represents  an 
ordinary  boundary  layer.  If  p(l)  >0  we  do  analogously. 

ii)  If  p(-l)  =  0  and  p* (-1)  <  0,  we  have  to  add  to 
the  set  of  trial  functions  a  function  of  the  same  type  as  the 
other  tJk,  whose  jump  is  concentrated  near  -1.  If  p(l)  =0 
and  p' (1)  =0  we  do  the  same  at  x  =  +1. 

iii)  If  p  has  a  multiple  zero,  this  zero  generates  a 
point-spectrum  that  becomes  dense  on  the  whole  positive  real 
axis,  if  e  tends  to  zero,  cf.  [11] .  More  specifically,  if 

(a)  ^  0  and  p(x)  =  0((x  -  ot)k)  (x  -*■  a)  ,  then  the  dis¬ 
tance  between  two  subsequent  eigenvalues  is  of  the  order 
0(eY),y  :=  (k-l)/(k  +  l).  If  p  changes  sign  at  a  with  a 
non-negative  slope,  the  smallest  eigenvalue  generated  at  this 
point  is  of  exponentially  small  order  and  it  has  to  be  taken 
in  account  as  before.  If  p  changes  sign  with  non-positive 
slope,  or  if  p  does  not  change  sign,  the  smallest  eigen- 

v 

value  is  bounded  away  from  zero  by  a  distance  of  order  0(e), 

hence  we  loose  only  an  extra  factor  e  ^  in  going  from 

(6.28a)  to  (6.28b).  However,  at  a  point  a  where  p 

changes  sign  a  boundary  layer  function  has  to  be  constructed 

and  added  to  the  trial  space?  the  boundary  layer  is  of  width 
1/ (k+1) 

C.  • 


i  .1 

ja 
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iv)  The  whole  analysis  carries  over  to  problems 

T 

governed  by  equations  of  type  -eAu  +  £  Vu  =  0  on  a  bounded 
domain  in  several  dimensions,  cf.  [6],  [10]  and  [12]. 

v)  The  estimate  (6.19)  yields  a  better  estimate  of  the 
first  eigenvalue  of  L£  than  the  formulae  of  Friedman  and 
Ventcel  and  Freidlin,  cf.  [13].  It  looks  even  not  too  diffi¬ 
cult  to  obtain  an  approximation  with  a  relative  error  of 
order  0(e)  by  considering  the  approximation  of  its  eigen¬ 
function  in  E  in  somewhat  more  detail. 
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